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Abstract 

A relativistic kinetic Fokker-Planck equation that has been recently pro- 
posed in the physical literature is studied. It is shown that, in contrast to 
other existing relativistic models, the one considered in this paper is invariant 
under Lorentz transformations in the absence of friction. A similar property 
(invariance by Galilean transformations in the absence of friction) is verified 
in the non-relativistic case. In the first part of the paper some fundamental 
mathematical properties of the relativistic Fokker-Planck equation are estab- 
lished. In particular, it is proved that the model is compatible with the finite 
propagation speed of particles in relativity. In the second part of the paper, 
two non-linear relativistic mean-field models are introduced. One is obtained 
by coupling the relativistic Fokker-Planck equation to the Maxwell equations 
of electrodynamics, and is therefore of interest in plasma physics. The other 
mean-field model couples the Fokker-Planck dynamics to a relativistic scalar 
theory of gravity (the Nordstrom theory) and is therefore of interest in grav- 
itational physics. In both cases the existence of steady states for all possible 
prescribed values of the mass is established. In the gravitational case this 
result is better than for the corresponding non-relativistic model, the Vlasov- 
Poisson-Fokker-Planck system, for which existence of steady states is known 
only for small mass. 



1 Introduction 

Fokker-Planck equations provide a continuous description of stochastic particles dy- 
namics. The most basic example is Brownian's motion, the stochastic motion of a 
test particle immersed in a fluid in thermodynamical equilibrium. Provided the test 
particle is much heavier than the molecules of the fluid, it is possible to approxi- 
mate the microscopic forces acting on the test particle by two driving mechanisms: 
diffusion and friction. The kinetic equation that describes the evolution of the dis- 
tribution function / for the test particle is the linear Fokker-Planck (or Kramers) 
equation [31]: 

dtf+p- V,/ = V p • (ppf + aV p f) . (1) 

The distribution function / is a non- negative function of the variables (t, x,p), where 
(x,p) are the phase-space coordinates (position and momentum) and t > is the 
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time variable. We assume that the mass of the test particle is one. The positive 
constants /3, a are the friction and diffusion parameters, respectively. The stochastic 
differential equations for the trajectory of the test particle associated to (1) are given 
by the system (4) in Section 2. 

Fokkcr-Planck equations like (1), or variants thereof, have several applications 
in different fields of physics and engineering. In astrophysics, for example, they 
model the effect of interstellar nebulas in a galaxy [32] or even dark matter [30]. 
In plasma physics, Fokker-Planck equations take into account the effect of grazing 
close encounters among the ions (the heavy particles) and the electrons. 

A questionable feature of equation (1) is that the diffusion term aA p f in the right 
hand side operates with infinite velocity: if the particles are initially distributed 
in a compact region of space, i.e., the initial distribution f(0,x,p) is compactly 
supported in the variable x, there will be instantaneously a non-zero probability 
(i.e., / > 0) to find particles everywhere in space. This property is incompatible 
with the well-established physical law that prevents particles from moving faster 
than light. Recent works in the mathematical and physical literature put forward 
two possible ways to eliminate this undesirable feature. One consists in replacing 
the classical linear diffusive (Laplace) operator with a non-linear diffusion term, 
as in the so-called "relativistic" heat equation, see [1]. A mathematically simpler 
solution is to replace (1) with a model that is still linear and, at the same time, 
consistent with the relativistic mechanics of particles, where the property of finite 
propagation speed enters in a natural fashion. The purpose of the present article is 
to begin the mathematical study of one such relativistic linear models. 

The physical literature abounds of proposals for what should represent the cor- 
rect relativistic generalization of (1), see for instance [10, 16] (and [12, 17] for an 
overview and an historical background to the relativistic theory of Brownian mo- 
tions) . Thus the first problem to face is the choice of the relativistic Fokker-Planck 
equation to consider. In this paper we pick the following equation: 

dtf+P-Vxf = V p -(pfp + aDV p f), (2) 
where p is the relativistic velocity, 

- = E 

and D is the relativistic diffusion matrix given by 

= I+p®p 

The previous model coincides with one of the equations proposed in [16], namely [16, 
Eq. (47)] and it is the subject of a recent series of papers by Haba [20, 21, 22, 23]. 
In these references several generalizations of (2) are introduced, including models 
for massless particles, for particles with spin and models with more general friction 
terms 1 . 

In Section 2 we justify our choice for the relativistic Fokker-Planck equation (2) 
by showing that it maintains certain important physical properties satisfied by the 
non-relativistic model (1). In particular we will show that for /3 = 0, i.e., in the 

x We are grateful to Prof. Haba for pointing out his work to us. 
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absence of friction, equation (2) is Lorentz invariant. Similarly, equation (1) is in- 
variant by Galilean transformations when (3 = 0. (Note that in both the relativistic 
and non-relativistic case the friction term breaks the equivalence of inertial refer- 
ence systems.) In Section 3 we prove that the solutions of (2) enjoy some other 
physically and mathematically desirable properties, in particular that they behave 
consistently with the finite propagation speed of particles. 

For the applications in astrophysics (rcsp. plasma physics), it is necessary to add 
the interaction of the particles with the self-generated gravitational (resp. electric) 
field. In the non-relativistic case this leads to the non-linear Vlasov-Poisson-Fokker- 
Planck system: 

dtf+p- V x / - V X U ■ V p / = V p ■ (pf + Vp/) , (3a) 
A x U = Xp, p(t,x)= [ f(t,x,p)dp, (3b) 

where we set all physical constant equal to one and where A = 1 in the gravitational 
case, while A = — 1 in the plasma physics case. In the second part of the paper 
(Sections 4-5) we introduce the corresponding relativistic model. In the plasma 
physics case we couple the relativistic Fokkcr-Planck equation (2) to the Maxwell 
equations of electrodynamics. The resulting model is the Vlasov-Maxwell-Fokkcr- 
Planck system. Note that this model is different from the one considered in [2, 28, 
36], which uses the non-relativistic Fokker-Planck equation (1). In the gravitational 
case we couple the Fokker-Planck dynamics to a relativistic scalar theory of gravity, 
the Nordstrom theory, which has already been used as a toy model for Einstein's 
theory of general relativity, see [4, 5, 6, 33]. Unfortunately there are fundamental 
difficulties, briefly recalled at the beginning of Section 5, in formulating a Fokkcr- 
Planck theory in general relativity. In this paper we prefer to avoid this issue and 
consider instead a toy model, which we call the Vlasov-Nordstrom-Fokkcr-Planck 
system. 

Our main result for the Vlasov-Maxwell-Fokkcr-Planck and Vlasov-Nordstrom- 
Fokker-Planck systems (with an external confining potential) is the existence of 
steady states solutions for all possible values of the mass. We do so by variational 
techniques inspired by [13]. Note that in the gravitational case our result is better 
than for the Vlasov-Poisson- Fokkcr-Planck system (3)a=i, for which the existence 
of steady states is only known for a properly small mass [3]. The main advantage 
of the relativistic model compared to the non-relativistic one is that the energy of 
the Vlasov- Nordstrom- Fokker-Planck system is positive definite. 

2 Derivation of the relativistic Fokker-Planck model 

A common way to derive Fokker-Planck type equations is to start from a system of 
stochastic ordinary differential equations (SODEs). The Fokker-Planck equation is 
the partial differential equation satisfied by the law of the stochastic process solving 
the SODEs. For instance in the case of the kinetic Fokker-Planck equation (1) the 
relevant SODEs are given by 

x(t)=p(t), p(t) = -pp(t) + V2aB(t), (4) 

where B(t) is the standard Brownian motion in M. d , i.e., a centered Gaussian process 
with covariance (B(t), B(t')) = S(t — t'), see [11, 27, 31] for details. Following this 
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approach to derive a relativistic Fokker-Planck equation is problematic for at least 
two reasons. Firstly it is not so clear how to define a "standard" relativistic Brow- 
nian motion. Secondly, there are multiple ways to derive a Fokker-Planck equation 
from a system of SODEs, which lead to different partial differential equations for 
the law of the stochastic process. For instance, equation (1) is obtained from (4) 
using Ito's calculus, whereas a different Fokker-Planck equation would be obtained 
by using Stratonovich's calculus 2 . As a consequence of these "ambiguities" , there 
exist different models in the literature which are named "relativistic Fokker-Planck 
equation", see [12, 17] for a review. 

The purpose of this section is to justify our choice for the relativistic Fokker- 
Planck model which will be studied in the rest of the paper. In particular we 
will show that it is possible to "derive" a relativistic Fokker-Planck equation by 
merely demanding that certain physical properties of the non-relativistic model be 
maintained in the relativistic case. We shall not refer in any moment to the SODEs 
for the (relativistic) stochastic process, although it will be finally observed that our 
equation coincides with one of the models derived in [16, 17] by stochastic calculus 
methods. 

We are interested in the following two important properties of the non-relativistic 
Fokker-Planck equation (1): 

(NR1) In the absence of friction, i.e. when /3 = 0, (1) is Galilean invariant 3 . This 
means that under the change of variables 

t = t , x = x — ut, p = p — u, fit, x,p) = f(t, x,p), 

f is a solution of (l),s=o if and only if / is a solution, Vu € R d . 
(NR2) The Maxwcllian distribution function 

M(p) = e -^l 2 /2- 

is a static solution of (1). In fact, up to a multiplicative constant, it is the 
only global equilibrium of the equation. 

We propose now a relativistic generalization of (1) by requiring that the relativis- 
tic analogues of the properties (NR1) and (NR2) hold for the new model. Precisely 
we require that the relativistic Fokker-Planck equation should satisfy: 

(Rl) Invariance under Lorentz transformations in the absence of friction, i.e., under 
the change of variables 4 



uo = v 1 + \u , t = uot —u - x , x = x — ut H : — 7T—u{u ■ x), 

m 

p = p-u^fl + \p\ 2 + U ° u(u-p), f{t 7 x 7 p) = f(t,x,p), 
M 

/ is a solution of the frictionless equation if and only if / is a solution, Vw S M d . 

2 These two difficulties are in some sense equivalent, since one can modify for instance the 
SODEs (4) to end up with (1) through Stratonovich's calculus. 

3 The friction term V p ■ (/9p/) breaks the Galilean invariance of (1), since it corresponds to the 
microscopic velocity-dependent force F = —f3p(t) in (4). 

4 We fix c = 1, where c is the speed of light. 
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(R2) The function g defined by 

dip) = e^^+^, 

must be a static solution, for some constant 7 > 0. 3 is known as the Juttner 
distribution (or relativistic Maxwcllian). 

The simplest and, in our opinion, most natural way to obtain (Rl) is the fol- 
lowing. Firstly we replace the transport term in the left hand side of (1) by its 
relativistic counterpart 

d 

VT+Wdt +P • = = 



with p° = ^/T+|pp, p = (p 1 , ■ ■ ■ ,p d ), do — dt and di = d x %. Secondly the diffusive 
operator A p = V p - V p on the right side of (1) is replaced by the Laplace-Beltrami op- 
erator A p h ^ over the Riemannian manifold (R rf , h), where h is the hyperbolic metric, 
i.e., the Riemannian metric induced by the Minkowski metric over the hyperboloid 
Sj = {(p°,p) ■ p° = \/i + \p\ 2 }- The fact that the operator A p is Lorentz invariant 
is clear, since the Lorentz transformation in the momentum variable corresponds to 
a translation over the hyperboloid $). The components of the metric h in the base 
dpi ® dpj of the linear space of second order covariant tensor fields on are given 
by 

hij &ij PiPj : 

where pk = Skip 1 and p = p/po is the relativistic velocity. Note that the position of 
the indexes (above or below) is changed using the Euclidean metric. Let (ft, -1 ) 4 - 7 = 
5** + p l p> denote the inverse matrix of h^, i.e., (h~ 1 y k hkj = 8j, and denote \h\ = 

det(hij) = (1 + |p| 2 ) -1 . The action of the Laplace-Beltrami operator Ap on scalar 
functions is given by 

AW/ = -±=d pi (VWKh-'fdpjf) . (5) 

Therefore the frictionlcss relativistic Fokker-Planck equation is 

a t f + p.V x f = ad pi ^^ dpj fj, (6) 

where a > is the diffusion constant. 

To achieve (R2) we add a friction term d p i(q l (p)f) to the right hand side of (6) 
such that the current 

A l =a JjLL dv f + q 1 f 

^TTW P 

vanishes for f = 3- It is straightforward to verify that this happens if and only 
if q l (p) = 707/, leading to the following relativistic Fokker-Planck equation with 
friction: 

dtf + P- V x f = Pfp 1 + o Z^ d v i f ), (7) 

V v 1 + \p\ J 



5 We will adopt the Einstein convention for the sum over repeated indexes. Greek indexes go 
from to d and Latin indexes from 1 to d. 



where j3 = is the friction parameter. 

Our purpose in the rest of the paper is to initiate the mathematical study of (7). 
Before proceeding, we modify (7) in two standard ways. Firstly, we set all physical 
constants to unity, i.e., /3 = a = 7 = 1; our results are independent from the value 
of the physical constants. Moreover, in order to guarantee the existence of finite 
mass equilibria in the whole space, we assume that the system is subject to the 
action 6 of an external confining potential V — V(x), and write the equation under 
study in the following final form 

Otf+P- y x f - y x V ■ V p f = V P (DV p f + pf), t > 0, p G R d , x G R d , (8) 

where D is the matrix D l: > = (S^ +p l p-')/ y/l + \p\ 2 (the diffusion matrix). Through- 
out the paper we assume V G C 1 and 

e^el'ft 3 ). (9) 

To conclude this section we remark that (7) coincides with one of the equations 
proposed in [16], namely [16, Eq. (47)]. In this reference the authors derive three 
different relativistic Fokker-Planck equations starting from a particular relativistic 
Langevin dynamics and using the pre-, mid- and post-point rule of discretization for 
stochastic integrals, see also [17]. Equation (7) is the only one, among the equations 
introduced in [16], that satisfies the properties (R1)-(R2) above. 

3 Basic properties of regular solutions 

In this section we prove some fundamental properties of regular solutions of (8). 
By regular solution we mean that 

< / G C([0,oo),L 1 (R d x R d )). 

Since the techniques we use are rather standard, some proofs will only be sketched. 

3.1 Cauchy problem 

We begin by sketching the proof of global existence and uniqueness to the initial 
value problem in the class of regular solutions. Let f- ln G L 1 denote the initial datum 
of /, i.e., f in (x,p) = f{0,x,p). 

Theorem 3.1. Given < f m G L , there exists a unique global regular solution. 

Proof. Approximate the external potential by a smooth function and the initial 
datum by a sequence fi n , m of smooth, non-negative functions with compact support. 
By the result proved in Appendix A, for each fixed m G N there exists a unique f m G 
C([0, 00), L 2 (diy)), solution of (8), where dv is the measure Av — cxp(-\/l + |p| 2 + 
V) d(p, x). Moreover by standard methods (see [9, 18, 35] for instance) one can prove 
the ^-contraction property: \\f k - fmWh 1 < ||/in,fc - /in.mlli 1 - Thus the sequence 

6 The action of the external potential is equivalent to that of a spatially dependent friction term, 
which can be seen by writing (8) in the form 

dtf + p-V x f = \7 p - (DV p f + f(v + V X V)). 
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f m converges in L 1 to a regular solution. The uniqueness is also a consequence of 
the /^-contraction property. The non-negativity of regular solutions can be proved 
by studying the evolution of a suitable regularization of sign(/) (see again [9, 18, 



We remark that it is possible to prove global existence and uniqueness of solu- 
tions with lower regularity, see [37] for the non-relativistic case. 

In the proof of the next results it will be assumed that the solution is smooth 
and decays rapidly at infinity. The generalization to regular solutions is achieved 
by introducing first a suitable smooth positive approximation f £ , for which the 
following calculations hold up to error terms that vanish in the limit toward a 
regular solution (i.e., e — > 0). We refer to [3] for the details of this procedure in the 
non-relativistic case. 

3.2 Finite propagation speed 

The first property that we want to emphasize is that equation (8) is compatible 
with the finite propagation speed of particles in relativity. 

Proposition 1. Assume that /; n = for \x — xq\ < to, where (to, Xo) € (0, oo) x M. d . 
Then / = for (t,x) € A(to,xo), where 



is the past light cone with vertex on (<o,-^o) oxid base on t = 0. In particular, if 



fin = for \x\ > R, for some R > 0, then f = for \x\ > R + t, for all t > 0. 
Proof. Introduce the density and the current density: 



Clearly \j\ < p and the continuity equation holds: dtp + V ■ j = 0. The result then 



35]). 



□ 



A(*o,xo) = {(*,a;)G [0,*o] x M rf : |a?-a | <*o-*} 




follows by Lemma B.l in Appendix B. 



□ 



3.3 Mass conservation and entropy identity 

Given a regular solution /, the mass is 




(10) 



and the free energy, or (relative) entropy functional is 



Q[f)(t)= [ f(t,x,p) ( VT+W +V(x) +\og f(t,x,p)) dpdx. (11) 
.raid \ / 



The next proposition studies the evolution of the functionals M, Q. 



Proposition 2. For a regular solution the following holds. 



(i) The mass is constant: M[f] = M[f- m ]. 
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(ii) If Q+[/in] < oo, where 

Q+[f]= f / ( V^+W +V(x)+ log+ /) dp dx, (12) 

Js 2d v ' 

log+/ = max(0,log/) ; then / log / e C([0, oo), L 1 (R 2d )), 

/ / ^ (p)0 p * ( y/JJs) dJjfjt) 3 dp dxds< oo 
and the entropy identity holds: 

\p)d p \^/JJd)d pl {yJU)3dpdx. (13) 



Proof. Proving the conservation of mass is straightforward. As to the entropy iden- 
tity (13), we begin by computing 



dQ 
dt 



I d t f(^T+W+V + \ogf)dpdx. 

jR 2d V ' 

We define d t f = FP[f] - T[f), T = p-V x - V ' X V ■ V p . First we see that 

(p ■ Vx/ - y/T+W^xV • V p / dp] dx 



T{f]^TTWdpdx = 

i 

= [ V x -{pf)dpdx+ [ p-V x Vfdpdx 

jR 2d Js. 2d 

= [ p-V x Vfdpdx. (14) 

Js? d 

For the integral of T[f]V we have 

T[f]Vdpdx= f (Vp-V x f -VW X V -V p f)dpdx 

i JK 2d 

= - f V x V-pfdpdx- [ V p -(VV x Vf)dpdx 

p-W x Vfdpdx. (15) 



For the integral of T[f] log / we use that for z = (x,p) and A a vector field such 
that V z ■ A = 0, there holds 

V,-[A(/log/-/)]=Alog/-V z / 

and therefore, taking A = (p, — W X V), we get 

f T[f](logf)dpdx= [ (p -X7 x f (log f) -V X V -V p f (log f))dpdx 

JK. 2d JR 2d 

= f Alogf-VJdpdx 

jR 2d 

= [ V z -[A(flogf- f)]dpdx = 0. (16) 

jR 2d 



Adding (14)-(16), we see that 

/ T[f] (^T+W+V + \ogf) dpdx = 0. (17) 
Now, for the term FP[-] we integrate by parts and obtain 
/ FP[f] i^T+W +V + log f) dp dx 

= I (Vl + W + V + log /) V p ■ (fp + DV p f) dp dx 
Pi + jdjf) {fp 1 + l)' l <)„ f) dpdx 



Using 



and 



we obtain 



/ 7 (fPi + d p , f) {fp 1 + D ij d pl f) dp dx. 

JVL 2d J 



rtj . 5v Pi +pi\p\ 2 j 
= 1 + =P 



^ (a/771) = — ^==(r7fc + r% k f), 



'[/] (Vl + W 2 + V + log /) dpdx 

z^'GO^ (y/71) d pl (yjjs) d dpdx. (is) 



FP[, 
= -4 

Adding (17) and (18) concludes the proof. □ 
3.4 Steady states 

Recall that e~ v <E L 1 . It is clear that, for each M > 0, there exists a unique regular 7 
static solution with mass M of (8), which is given by 

M 

fo(x,p) = m M {x,p) = -q3v{x,p), (19a) 
where dv{x,p) = e ~(^ 1+|p|2+v ) and 9=7 3 v (x,p) dpdx. (19b) 



Moreover, as in the non-relativistic case, one can prove that the equilibrium 
solution is a minimizer of the entropy functional. To sec this, we first recall the 
following general result proved in [13, Lemma 1.1], which will also play a crucial 
role in the following sections. 



7 In Appendix A it is proved that the operator L, defined by writing the equation (8) in the form 
dtf = Lf, is hypoelliptic, provided V 6 C°°. From this property one obtains that the equilibria 
of (8), which solve Lf = 0, arc automatically smooth. 
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Lemma 3.2 ([13]). Let us consider f2 C M. d measurable and the functional 

H[g] = I g(y) log g(y) dy + f g(y)h(y) dy, (20) 



with g G L 1 (f2) non-negative such that g(\ogg) + G L 1 ^). If h G L (Q; g(y)dy) is 
such that e~ h G L 1 (f2; dy), then g log g £ L 1 (Cl; dy) and 



H[g}-H[m g ] > 



where rn g (y) 



Jnddy u 



Ia e ~ h d y 



An immediate consequence of the previous lemma is a characterization of the 
minimum of H: 

Corollary 1. With the same hypotheses of Lemma 3.2, 

H(M) = inf ^U[g] : g > 0, g G L 1 ^), J^g(y) dy = M 

is bounded from below for any M > and 

M 



H{M)=n\g] = M log 



J n e h dy 



with g = M 



J Q e h dy 



In fact, g is the only minimum ofH(g). 



If we take N = 2d, y = (x,p), Ct = R 2d , g = f and h = y/l + \p\ 2 + V, we have 
H(g) = Q(f) and g = n%M- Thus we obtain 

Corollary 2. Assume that f G L 1 (R 2d ), f > are such that 

Q+[f}= [ f( v/i + bl 2 + V(x) + log+ A dp dx < oo 

and e~ v G L 1 (R d ). Then, 

Q[f] - Q[m M ] > \ f M (Vf&P) ~ V™m(x, P )Y dpdx 
and niM is the unique minimum of 
Q(M) = inf f Q[/] : / > 0, / G L 1 (M. 2d ), f f(x,p) dpdx = M, Q+[f] < oo 

= Mlog 



M 



J R2d e~(V^W+v)dpdx 



In the next sections we shall generalize this result to the non-linear Vlasov- 
Maxwell-Fokker-Planck and Vlasov-Nordstrom-Fokkcr-Planck systems. 
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4 The Vlasov-Maxwell-Fokker-Planck system 



In the present and next sections we consider two non-linear mean field models built 
on the relativistic Fokker-Planck equation (8). These models provide a relativistic 
generalization of the Vlasov-Poisson-Fokkcr-Planck system in the plasma physics 
case (present section) and in the gravitational case (next section). For simplicity 
we shall consider only the three dimensional case, i.e., x,p £ M. 3 (the field equations 
change with the dimension). 

The relativistic model for plasmas is obtained by coupling the relativistic Fokker- 
Planck equation 

dtf + P ■ V x / + F ■ V p f = d p , (f P l + D ij d p j /) , (21a) 
for the Lorentz force field (with external potential) 

F : [0, oo) x R 3 x R 3 — > K 3 , F = E+pxB-V x V (21b) 
and the system of Maxwell equations given by 8 



with 



d t E = V x AB-j, (i) V x -E = p, (u) 
d t B = -V x AE, {til) V x -B = 0, (iv) 

p(t,x) = / f(t,x,p)dp, j = I pf(t,x,p)dp, 



(21c) 



(21d) 



where E, B : [0, oo) x I 3 -> I 3 are functions of (t, x) that represent the electric and 
the magnetic field, respectively. Note that (p, j) satisfies the local conservation of 
charge 

dtp + V, ■ j = 0, (22) 

as a direct consequence of (21a), which makes it consistent to couple the Maxwell 
equations and the Fokker-Planck equation. 

The system (21) will be called the (relativistic) Vlasov-Maxwell-Fokker-Planck 
system, or VMFP for short. It generalizes the Vlasov-Poisson-Fokker-Plank (VPFP) 
system in the plasma physics case. Therefore (21) takes into account relativistic 
effects in a plasma, such as the propagation of electromagnetic waves. We remark 
that there exist other models in the literature which are named "Vlasov-Maxwell- 
Fokker-Planck" , see [2, 28, 36]. These systems couple Maxwell's equations to the 
non-relativistic Fokker-Planck equation (1). 

This section continues by proving the mass conservation and the entropy iden- 
tity of time-dependent solutions and the existence of steady states to VMFP. The 
analysis of time-dependent solutions is only formal, since there is no proof of the 
existence of solutions with enough regularity to which apply the argument below. 
We shall use the terminology "regular solution" of VMFP in a loose sense, meaning 
that the solution is non-negative and sufficiently regular to enable the following 
calculations. 



'Up to a suitable normalization of the physical constants. 
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4.1 Formal properties of regular solutions 

The mass of regular solutions of (21) is defined by (10); the entropy functional is 
defined as 

K\f,E t B] = [ f(y/TTW+V + logf) dpdx + ± [ (\E\ 2 + \B\ 2 )dx. (23) 

JR* k J I J R 3 \ J 

Proposition 3. For regular solutions of (21) we have: 

(i) The mass is preserved: M(t) = const. 

(ii) The entropy functional satisfies 

dJC _ 
~dt ~~ 



(p)d p <(^f7d)d pl (^J/d)dd P dx. (24) 



(Hi) Let e v E L 1 (M 3 ). Regular static solutions of (21) with mass M verify 

(f (x,p),E (x),B (x)) = (m M (x,p),-VU(x),0), (25a) 

where 

m M (x,p) = ^r3v(x, P )e- uix \ G= f e- u ^3v(x,p)dpdx, (25b) 
H Jr® 

and U is a solution of 

— AU = p, p = / m M (x,p) dp. (25c) 
Jr 3 

Proof. Proving (i) is straightforward. To achieve (ii) we write K, = Q+I, where Q 
is given by (11). Thus 



dJC dQ dl f / r \ dl 

* = ~t + Tt = J tf [S^W + v + lo ^J d ^ d * + Tf 



Let d t f = FP[f]-T[f], where in this case F = E+px B -V X V for T[]. Therefore 
we only need to calculate the derivative of I[E, B) and the part of dQ/dt containing 
the term E + p x B in T[-], since the other terms from T[-] and FP[-] are the same 
as in the linear case, cf. Proposition 2. Using (21ci) y (21ciii), we have 

^ = f (E-d t E + B- d t B) dx 
dt J M 3 

= J (E ■ (V, A B - j) + B ■ (-V x A E)) dx 

= f [ (E ■ (Vx A B) - B ■ (V x A E)) - E ■ j] dx 

Jr 3 

= [ V x ((BxE)-E-j)dx = - [ E-jdx. 
Jr 3 Jr 3 



Moreover 



/ (E + pxB)-V p f(\ogf + V)dpdx=[ V p - [(E + px B)(f\ogf-f)]dpdx 

JR e il« 

+ / V p • (V(E+p x B)f)dpdx = 0. 

JR 3 
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We also have 

/ V 1 + \p\ 2 (E+P x B) ■ Vpfdpdx = - / p- (E+px B)fdpdx 

E ■ j dx. 

The second equality is due to the orthogonality between p and p x B and the 
definition of j. The proof of (ii) follows easily. As to (hi), we hrst notice that from 
(21ciii) and since dtBo = 0, we have = dtBo = — V A Eq =>■ 3 U{x) such that 
Eq = — Vt/(x). Using (21cii) we obtain —AC/ = p. Moreover by (24) applied to 
static solutions we observe that fo(x,p) = a(x)3(p) for some non-negative function 
a = a(x). In particular, j = (since it is the integral of an odd function) and the 
equations for the field Bo are equivalent to V x Bq = V • Bq = => £?o = 0. Now 
replacing / = &d, E = — VC/ and B = in (21a) we obtain 

Va + ap- V(U + V) =0. 

It is clear that the only non-trivial regular solution of the previous equation is 
a = Ce~ u ~ v , where C is any positive constant. The value C = M/d follows by 
the definition of M. □ 




4.2 Existence of steady states 

In this section we prove the existence of (regular) static solutions for the system 
(21). In particular, we want to show that the free energy functional 

JC[f,E,B] = Q[f}+T[E,B) 

= f f(V^+W+V + \ogf)dpdx + ±[ (\E\ 2 + \B\ 2 )dx, 

JR6 V / Z J R 3 \ / 

subject to 

V-E = p, W-B = 0, j fdpdx = M, 

attains its minimum exactly in the static solution of (21) with mass M. The fol- 
lowing proof generalizes the one given in [13, Prop. 2.2] for the VPFP system. Note 
that the variational problem for VMFP differs from that of VPFP studied in [13] in 
two aspects. Firstly, the electromagnetic field appears as an independent variable in 
the entropy functional, while for VPFP the electric field is given by the convolution 
product of p with l/(47r|a;|). Secondly, in the variational problem for VMFP there 
appear the local constraints V ■ E = p, V ■ B = 0. Nevertheless we will be able to 
reduce the problem at hand to the equivalent one for the VPFP system considered 
in [13]. In particular we will show that the above minimization problem is equiva- 
lent to minimizing a reduced entropy functional /C re( j that resembles the free energy 
in the non-relativistic case. To this purpose we use the following simple result. 

Lemma 4.1. The solutions of the variational problem 

inf 11(h) = inf / \h\ 2 dx, 
he's hes J R 3 

where D = {h e L 2 (R 3 ) : Vh = g}, g € L^IR 3 ), are of the form h = -VU , where 
-AC/ = g. 
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Proof. Let be a test function. The first variation of 1Z evaluated on a critical 
point has to vanish, which implies 

d f d f 

—K(h + tcj>) |t= = / -t: \h + t<j)\ 2 | t=0 dx = / 2fc-^dx = 0. 

at J R 3 at J R 3 

In particular, we can consider test functions of the form (f> = V A v, which entails 
= / h-V Avdx = - V A h ■ v dx, 



for all v <G C^°(M 3 ). From here we infer that V A h = and as a consequence, there 
exists U such that h = —VU. Substituting this value in V • h = g concludes the 
proof. □ 



Next we define 



J f (Vi + bl 2 + \u + v + log/) dpdz, 

with — AU = p and p = / R3 / dp. 

Proposition 4. Recall the definition (12) of Q+[f]. Let 

K(M)=ivf[lC[f,E,B] : / > 0, / G L 1 (R 3 ), ||/|| Ll(RM) = M, Q+[/] < oo, 

(£7, S) e L 2 (M 3 ) x L 2 (R 3 ) , V • E = p, V • B = o| 

and assume that e~ v € L 1 (M 1 ). Then, 
(i) K(M) = inf {/C red (/) : / > 0, / e L 1 ^ 6 ). 11/ IIl^) = M > C+t/1 < 
(zi^ K(M) is bounded from below for any M > 0; 

fmj JTie minimizer is unique and is given by (25), i.e., K(M) = /C re d (tim)- 
Proof. To show (i), let X denote the minimizing space and define 

*i = {/ e ^(R 6 ) : / > 0, ||/|| L1(R6) = M, Q+[/] < oo} , 
X 2 = {(E,B) e L 2 (R 3 ) x L 2 (R 3 ) : V • £ = p, V • B = 0} . 
The minimum (if it exists) verifies: 

K(M) = inf {£(/, £, B)} = inf {inf {J(£, B)} + Q(/)| 

= inf(l / IVtffds + GC/)} 



=ff dp 



inf {/C red (/)} 



since by Lemma 4.1, for g\ — p and 32 = 0, we have E = — Vf7, — AJ7 = p and 
_B = —WU, —AU = 0, which implies {7 = 0. On the other hand, we see that 

- / \WU\ 2 dx=- I -UAUdx=- I pUdx = j -fUdpdx 

2 JR3 2 J R 3 2 J R 3 J R 6 2 
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and the original problem is therefore reduced to minimize the functional /C re d(/), 
which, up to substituting yl + \p\ 2 with |p| 2 /2, coincides with the free energy in 
the non-relativistic case. Thus the claims (ii) and (iii) can be established as in [13, 
Prop. 2.2]. □ 

To conclude this section we remark that the existence of steady states to the 
VMFP system can be established also by studying directly the equation (25c), as 
done in [19] for the non-relativistic case. The non-existence results proved there 
when e~ v ^ L 1 (see also [13]) are valid in the relativistic case as well. 



5 The Vlasov-Nordstrom-Fokker Planck system 

In this section we introduce yet another new model, which represents a relativistic 
generalization of the VPFP system in the gravitational case. It would be desirable 
to obtain such a model in the framework of general relativity, since the latter is the 
physically correct relativistic theory of gravity (as far as we know...), but this would 
lead inevitably to face fundamental difficulties. In fact the consistent modeling of 
dissipativc systems in general relativity is not yet understood, not even at a formal 
level, the main reason being that the Einstein equations by themselves imply that 
the mass/encrgy/momentum of the system must be conserved 9 . To overcome this 
(still unresolved) fundamental issue, instead of general relativity we shall use an 
alternative relativistic theory of gravity, the Nordstrom theory, which has already 
been used in the collisionlcss case as a toy model for the more complicated Einstein- 
Vlasov system [5, 33]. The resulting system — the Vlasov-Nordstrom-Fokker-Planck 
system — will be derived using an argument similar to the one applied in Section 2. 



5.1 Derivation of the model 

While for the VMFP system the background space-time is given by the manifold 
(M. 4 ,r]), where r\ is the Minkowski metric, in the present case we assume that the 
space-time is given by the Lorentzian manifold (U. 4 ,g), where 

9 = e Vr n, 

with (j> : M 4 — >• R a scalar field, which will play the role of the gravitational field. 
Let (t ) be a system of coordinates which set the Minkoski metric in the 

canonical form f\^ v = diag(— 1, 1,1,1). Then 

g = -e^dt 2 + ( -'>), , d.r'd.r '. (26) 

The geodesies of the metric (26) are the solutions of the following system of 
ODEs: 

dt n dx l i dp 11 „ „ „ . . 



9 The situation is similar to what happens in electrodynamics, where the Maxwell equations 
alone imply the conservation of charge (22) and therefore the dynamics of the coupled matter model 
must be compatible with it (which is true for the relativistic Fokkcr-Planck equation considered 
in the previous section). 
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where s is the geodesic parameter and T^ a are the Christoffel symbols of g: 

Ka = (dug<n + fan - diQuo) = Sid a 4> + 5%d v cj> - e~ 2 * 'g vaV T 'd y <f> ■ (28) 

Let us consider a system of particles with unit mass that move along the geodesic 
curves. The geodesic motion reflects the physical property that the particles inter- 
act only through the gravitational field. If we want to interpret p^ as the four- 
momentum of the particles, we need to impose that p M has length equal to —1, i.e., 
g^uP^p" = — 1. This entails 

p° = y^ + H 2 , \ P \ 2 = 6 ijP y. (29) 

Let f(t,x,p), x = (x ,x , x ) and p = (p ,p ,p ), be the distribution function 
of particles in the position x at time t and with four- momentum p^ = (p°,p) = 
(y 1 e~ 2 ^ + \p\ 2 ,p). Having assumed that the solutions of (27) are the particles tra- 
jectories, we obtain that / satisfies the equation 

P°d t f+ P -v x f-ri 1/P ^d pi f = o, 

where p° is given by (29). Substituting (28) in the last equation we obtain 

p°d t f + p ■ W x f - [2 (podtcj) + p ■ V x <t>) p + e-^V^] ■ V p / = 0. 

The previous equation is the Vlasov equation for collisionless particles. For the 
Fokker-Planck equation we need to add a diffusion and a friction term in the right 
hand side. Motivated by the discussion in Section 2, for the diffusion term we 
pick Ap 1 -*/, where h is the metric induced by (26) over the hyperboloid p° — 
\J e~ 2 ^ + \p\ 2 . It can be verified that 10 

hj = e 2 * (sy- e _ 2 ^| p | 2 j , Pi = hP j - 

We have 

\h\ =det/i= (e-^ + bl 2 )"^ 4 ^, (h^ 1 )^ =e- a *<y tf +pV. 



Therefore, 



- \ P \ 2 d pi I - - " ' pp d vi f . 

We then obtain the Fokker-Planck equation in the absence of friction in the following 
form: 



Sf 



2S<f>p- 



0,-20 + | p |2 



V p / = d pi = -d pj f , 30a 

V V e * + p J 



1() Although the metric g is not Euclidean, we keep using the metric 5ij for moving up and down 
indexes. 
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where 



Su = d t u + - = ■ V x u. (30b) 

Ve- 2 * + \p\ 2 



For the scalar gravitational held <f> we postulate the non-linear wave equation 



□0 := dfcf, - A x cj) = -e^ f f{t ; X ' P) dp, (30c) 

Jr3 Ve- 2 ^ + \p\ 2 

which has been justified in [4]. Now, doing the change of variables f(t,x,p) = 
f(t,x,e~ 2 ^p), the system (30) takes the form 

dti + v P + \ P \^ . v,/ - yj* + \ P \^ ■ vj = o pi (^( P )d p j) , 

(31a) 

Ucp = -e 2 * [ ~ f[t ' X ' p) dp, (31b) 
7m3 Ve 20 + \p\ 2 



where 



e ±<t>§ij + e^pipJ 



The system (31) is the Vlasov-Nordstrom-Fokker-Planck system in the absence of 
friction. It is invariant under the Lorentz type transformations given in [7]. To 
introduce a friction term, we first notice that for any given time independent scalar 
function 0o = 4>o(x), the left hand side of (31a) vanishes for 



f = fo(x,p)=e-V^°+\p\ 2 . (32) 

This suggests to introduce a friction term of the form V p • (qf) on the right side 
of (31a) such that 



Aj'(p)fl p< /+«7=0, if / = e-\/« 2 *+H a . 

It can be verified that q = e 2 *p. Adding this friction term and an external potential 
to (31a), we get 

dtf + V p \^fi?* + \p~\ 2 ) ■ V,/ - W x ^e^ + H 2 + V(x)j ■ V p / (33a) 

= d pi (Af(p)^/ + e 2 V/), 

Udp = -e 2 * I f{t ' X ' P) dp, (33b) 
Jr3 ^e 2( t> + \p\ 2 

where A 7 (p) is given by (31c) and where we removed the tilde for notational simplic- 
ity. The system (33) will be called the Vlasov-Nordstrom-Fokker-Planck (VNFP) 
system. 
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5.2 Formal properties of regular solutions 

Besides regularity, we assume that the solutions of VNFP are such that is 
bounded in any finite interval of time. This is true as soon as the initial data 
for the field equation (33b) are bounded. To see this, note that regular solutions 
of (33b) verify <fi = </>hom + where 0hom solves the wave equation dehorn = 
with the same data of (f> and ip solves (33b) with zero data. Since the right hand 
side of (33b) is non- negative, then ip < 0, and therefore = e^ hom e^ < e^ hom is 
bounded, as we claimed. 

The mass of regular solutions of VNFP is defined by (10). The entropy functional 

is 

JC[f, <t>, d t 4>] =Jf (je* + b| 2 + V(x) + log f) dpdx 



= Q[f,4>]+l[4>,d t 4>]. (34) 

Note that the energy part of the entropy functional is positive definite, in contrast 
to the case of the gravitational VPFP system. 

Proposition 5. For regular solutions of (33), we have: 

(i) M(t) = constant. 

(ii) The entropy functional satisfies 

dJC 



= -4 
dt 



A; j (p)fy ( yfw*) d P , ( yfw*) t dp dx, (35) 



where ^{x,p) = e -V e2 *+H 2 . 
(Hi) Let e~ v G L 1 . Static solutions of VNFP with mass AI > are of the form 

(fo(x,p),<j>o{x)) = (m M (x,p),(l>o(x)), (36a) 



m M ( 



r(x,p) = %3v(x,p), dv = e-V^°+\p\ 2 -V i Q = f dv ( x , p ) dpdx, 

(36b) 

and (f>o solves 

A(f>o=e 2, f ™ M {X, P ) ^ 

JW3 V e20 ° + \P\ 2 

Proof. The proof of (i) is straightforward. To show (ii) , we first observe that 

^ = / d t f (Je** + W + V + \og f) dp dx + f ff^f.. dp dx. (37) 



11 Note that 3v S L 1 ^ 6 ), because is bounded. 
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Again we split d t f = FP[f] - T[f], where T = V x {y / e 2 't> + \p\ 2 - V) ■ V p 



V p (i/e 2 ^ + \p\ 2 ) • V-e. For the integral containing T we have 

/ T[f]Je*<t> + \p\*dpdx = - [ fj*^ dpdx, 
Jr» v Jr<> ^/e 2, t> + \p\ 2 

[ T[f]\ogfdpdx = Q, 

/ T[f]Vdpdx= [ p ' VxV dpdx. 

jk 6 jr8 v e2 * + bl 2 

Thus the term involving T[f] gives no contribution. Moreover 
J FP[f] (y/e 2 <t> + \p\ 2 + V + log /J dp dx 

= ^ (pi + jd pi f) (e^fp* + A«(p)^/) dpdx 

= -4j swjtodpi ym^d pj yjjd^j dpdx, 

where we used that 

e 2 V = A«(p)ft and d pk (y[f/9*} = ^= (f Pk + d pk f) 

-77=/ (d t ct>d?ct> + V x (t>-Vd t 4>)dx= [ dt^dx, 



On the other hand, 



which cancels the last term in (37) due to the field equation (33b). This concludes 
the proof of (ii). For the last statement, we use that, by (ii), static solutions must 
have the form /o( x, p) = a(x)d^° (x,p). Substituting in (33a) we obtain the equation 
p ■ (Vq + aSJ x V) = and therefore a = Ce~ v . □ 

5.3 Existence of steady states 

The existence of steady states for the VPFP system in the gravitational case is 
not yet well-understood. We mention that a small mass result is proved in [3] 
for the VPFP system using a fixed point argument inspired by [14, 15]. This 
argument applies mutatis mutandis to the VNFP system: Consider the equation 
for the gravitational potential of steady states, eq. (36c), which we rewrite in terms 
of u = —0o as 



e~ v Me" 2u f e"V e 2 "+H 2 



f Ke e-Ve- 2u +W 2 -V dpdx V y/e- 2u + \p\ 

Define by K the solution operator of (38), i.e., the convolution of the r.h.s. with 
l/(47r|x|). By standard estimates one can prove that, for M small enough, the 
operator K is a contraction in the space X = {v € L°°(R 3 ) : < v < 1} and so by 
the fixed point theorem we have the following result. 
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Proposition 6. There exists Mq > such that, for all M < Mq, the equation (38), 
with the boundary condition lim|a.i_ > . 00 u = 0, has a unique solution u G L°° . This 
solution defines, through (36), a steady state of the VNFP system. 

However in the relativistic case we can do much better, and prove existence of 
steady states for all masses. Let us denote 

r A/ = {/:R 6 ^M : / G L 1 ^ 6 ) , \\f\\ L i m =M, Q+[f] < co}, 

where Q + is dchned by (12), and recall that the space D 1 (R 3 ) is dchned as 

D 1 (R 3 ) = {4> G iio^K 3 ) : V</> <G L 2 and 4> vanishes at infinity}, 

where the condition of </> vanishing at infinity means that the set {x G M 3 : \4>{x)\ > 
a} has finite (Lebesgue) measure, for all a > 0. Functions in the space Z? 1 (M 3 ) 
satisfy the Sobolev inequality 

H\\l><v\\V<I>\\l', ^=^f~ V ^ (39) 

see [29, Thm. 8.3]. 

Theorem 5.1. For all M > there exists at least one solution 12 cf>o of (36c). 
Moreover the corresponding steady state, given by (36b), is a minimizer of the 
entropy functional: 

K(M) = inf{/C(/, 0, i>) , f G T M , <t> G D^R 3 ) , V G L 2 {R 3 )}, 

where K, is defined by (34), i.e., K(M) = K.(mM,4>o,0). 

Proof. First we notice that 

K(M)= inf £{f,4>), 

where £(f, <fi) = K{f, <f>, 0). We divide the proof in five steps. 
Step 1: K{M) is bounded. We have 

S(fA)> I f(\p\ + V(x)+logf)dpdx. (40) 
Jr 6 

Using Lemma 3.2 with g = /, h = \p\ + V, Q = R 6 we get 

Step 2: Weak convergence of minimizing sequences. Let (/ n ,</>n) be a minimizing 
sequence. Since </>„ is uniformly bounded in D 1 , and by the Sobolev inequality (39), 
there exists a subsequence, still denoted by </>„, and 4>q G D 1 such that 

4> n ^4>o in L e and V x 4> n -± V^o in L 2 . (41) 

12 By Proposition 6, the solution is unique for M small. 
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Next we establish the weak convergence of f n in L 1 by using the argument in [13, 
pag. 129]. Let us show first that /„ does not concentrate. If it did, we could find 
e > 0, a bounded sequence x n G R 3 and a sequence R„ — > oo such that 

/ fn(x,p) dpdx = e, for all n G N. 

J\x„-x\<R„ 

From (40) and Lemma 3.2 we have 

S(U,<M > [ /n(log /„ + \p\ + V(x)) dpdx 

J \x-x n \>R n 

f n (\ogf n + \p\+V{x))dpdx 

x n —x\<.R n 



>(M-e)log 



M-e 



f e-\p\~ v dpdx 

J \x„ — x\>R n 1 j 



H ^-,<J-'"'- v ^ )- <42) 



Since e ' p ' ''gL 1 , we have 

lim / e-^- v dpdx = 0, 



and lim / e _|p| - v dpda; = He" 1 "'-^!^™) 

and so (42) implies £(f n ,<f>n) —> oo as n oo, which contradicts the fact that 
{fn,4>n) is a minimizing sequence. Now we prove that /„ is tight. If not, we can 
find £ > such that, for all Rq > 0, there exists R > Ro such that 



lim / /„ dpdx > £. 

n ^°° J\x\ + \p\>R 

Whence, using again (40) and Lemma 3.2, 



£(fn,4>n) > / fndpdx 
\J\x\ + \p\>R / 



log / fndpdx 

J\x\ + \p\>R 

-log f e- M ' v dpdx 

J\x\ + \p\>R 



and so 

lim / e-^- v dpdx > £e~ K{M)/e > 0, 

R ^°° J\x\ + \p\>R 

which contradicts the fact that e~' p '~ v € L 1 . We conclude that there exists /o G L 1 
and a subsequence /„ such that 

U^foinL 1 . (43) 

S'iep 3: Pointwise convergence of minimizing sequences. As proved in [29, Cor. 8.7], 
the weak convergence (41) implies that 

4>n — > 0o, pointwise a.e. (44) 



21 



again up to the extraction of a subsequence. Moreover, by the argument used in [6, 
Lemma 5], we may assume that (f> n < almost everywhere. Next we show that 
f n converges pointwise a.e. (up to subsequences). Given a minimizing sequence 
(/„, 4> n ), define 



M 



3 - v / e 2 *"+|p| 2 -V 



f R6 e-V^W-v dpdx 

By Lemma 3.2 we have 

£{fm <t>n) - £{m n ,4>n) > ^ / ( \//n - \f™hiY dp dx. 



This implies from one hand that (m„,</>„) is again a minimizing sequence and, on 
the other hand, that lim n _ >00 (/„ — m n ) = pointwise a.e. after extracting a suitable 
subsequence. Moreover, since 

e -V^+\p\ 2 -v _> e-V^W-v^ pointwise a . e . 
and (by dominated convergence) 

e -V^"+|p| 2 -v dpda ,^ I e -Ve 2 *o+| P ! 2 -v dpda , ; pointwise a.e. 



f n — > ; e - x /e 2 *o+| p |2-y pointwise a.e. 

In particular we notice that Jo is strictly positive and bounded. 

Step Jf.: (/o, <po) is a minimizer. We prove that £ is weakly lower semicontinuous. 
Clearly 



liminf / \V x <j) n Y&x > / V^o \ dx. 

rwoo J J 

Moreover, by Fatou's lemma, 

liminf f f n (Je 2 ^ + \p\ 2 + V + log /„) dp dx > f f {Je^° + |p| 2 + V + log/ ), 

and the claim follows: K(M) = £(fo,4>o). 

Step 5: (/o, 0o) * s a steady state of the VNFP system. Since we already proved 
in step 3 that 

fo 



A/ ^/ e 2* 0+ | p |2_y 



J R6 e-V^ +bl 2 -^dpdx 

we only need to show that </>o solves the non-linear elliptic equation (36c). To this 
purpose we define <ph — 4>a + hi], where rj = rj(x) is any C°° function with compact 
support and h £ R. Using that < /o < oo and (j>o < 0, it is straightforward to 
show that £(fo,(j)h) is differentiable in h. The derivative at h = must vanish and 
this entails that 4>q solves 

A0 o = e 2 ^ / ; f ° dpdx 
in the sense of distributions. This completes the proof of the theorem. □ 
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A The Cauchy problem for the Fokker- Planck equa- 
tion on a Riemannian manifold 

In this appendix we discuss the initial value problem for the equation (8). In fact we 
shall study the problem for a more general equation than (8), where we allow for a 
general (positive definite) diffusion matrix D, a general velocity field (with non-zero 
gradient) in the transport operator and a general friction potential. Precisely we 
shall consider the initial value problem for the following equation: 

d t h + v(p)- V x h-V x V- V p h = A^h + Wh, t > 0, x e R d , p £ R d , (45) 

where Ap 9 ^ denotes the Laplace-Beltrami operator of a Riemannian metric g on 
see (5), and W,v are the vector fields 

Wh = g~ l V p \ogu ■ V p h, u = ^/deTg' e~ E , v = V p E, (46) 

for some non-negative function E = E(p). Equation (8) can be written in the 
form (45) by setting g = E(p) = yj\ + \p\ 2 , f = er E ~ v h. We prove the 

following 

Theorem A.l. Assume g,E,V € C°° satisfy e~ E , e~ v € i 1 ^ 3 ) and, for all 
p e K d , 

dct(d p iVj) 7^ 0, 

dpiig^dpjE) < ui, for some ui > 0, 

g ij d p iEd p3 E > 9\V p E\ 2 , for some 9 > 0, (47) 
where g lJ is the inverse matrix of g, i.e., g lJ gjk = S l k . Furthermore we assume that 
a ij (v) 

^pp^O, as \p\^oo Vi,j' = l,...d. (48) 

Then for all h- m e C c 1 (M <i x R d ) there exists a unique 

he C([0,oo),L 2 (d/i)), 

solution of (45) with initial datum h m , where /i is the probability measure defined 
by 

dfj, = Q- 1 e- E - v dpdx, 6 = / e-^-^dpdx. 

JR 2d 

For the proof we need the following lemma. 
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Lemma A. 2. Let A = -L + T, where L = A p g) + W andT = v(p) ■ V x - V X V- V p) 
and /,/i,/ 2 € C°°. T/ien 

(oj f md hLhe~ E dp = - j^g 13 d p ihd pj h e~ E dp; 

(b) f R2d hThdn = 0; 

fcj A{hf 2 ) = hAh + M/i - 2g*id^h d pl f 2 . 
Proof. The proof of (a) can be found in [8] . For the second statement we use that 

/ h(v(p) ■ W x h- W X V- W p h)e- E - V dpdx 

JR 2d 

= I (-V K /i- V p Eh + Vph- V x V)he- E - y dpdx, 

using that v(p) = T/ p E and integrating by parts. The proof of (c) follows by 
Leibnitz's rule. □ 

Proof of Theorem A.l. We generalize the proof of [24, Prop. 5.5], where the fol- 
lowing argument is applied to the non-relativistic Fokker-Planck equation (1) (with 
external potential), and the proof in [8, App. A], which studies the Cauchy problem 
for (45) when x € T d (the d— dimensional torus) without external potential. 
Denote TL = L 2 (d/i). Let us consider the operator 

A = v(p) ■ V x h - V X V ■ V p h - A p 9) h — Wh = T — L 

defined on D(A) = C c °°(]R 2d ). Equation (45) takes the form d t h + Ah = 0. Our goal 
is to show that the closure of the operator A generates a contraction semigroup on 
TL. To this purpose it suffices to prove that A is accretive and that the range of 
A + XL is dense in TL for some A > 0, see [24, Sec. 5.2]. 

That A is accretive follows by (a) and (b) of the previous lemma: 



(h \ Ah) n = -(h \ Lh) n + {h\Th) H = J g^dp.hdpjhdfi > 0. 

Next we show that A is hypoelliptic. Let a = \J <? _1 , the positive definite matrix 
such that a 2 = g~ l . A direct computation shows that 

d 

-A = Y,Y?+Y Q , 

i=i 

where Yq, Yi denote the vector fields 

Y Q h = (div p a) • aVp/i - g l] d pl E d p jh - Th, 
Yih = a^d p kh. 

In order to prove that A is hypoelliptic, it is enough to show that —A satisfies 
a rank 2 Hormander's condition, i.e., the vector fields Yj and Zi := [lo,Fi] form a 
basis of R 2d , see [25]. Observe that 

Zi = B^dpk + Cjd x j, 
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where G\ = a^dpkv^ and B is a d x d matrix whose exact form is irrelevant for what 
follows. Thus we can represent the linear transformation {d x i, d p j} — > {Yfc, Z{\ by 

/ a 
\ C B 

whose determinat is |detF| = dctaldctCI = dot g\dct(d p kv : ')\, which is positive 
because dct(d p iVj) is non-zero by assumption. Therefore, {Yj,, Zj} is a basis of R 2d . 

Finally, we prove that the range of A + A is dense in H for some A > 0. If h £ H, 
we must show that if 

{h\(X + A)f) n = 0, for all / G D(A), (49) 
then h = 0. Equation (49) is equivalent to h being a distributional solution of 

(A - L - T)h = 0. 

Since the operator in the left hand side of the latter equation is hypoelliptic, then we 
can assume h £ C°° . Now setting fa = </>, fa = <f>h in (c) of Lemma A. 2, multiplying 
by h, integrating and using that (h\(X + A)(4> 2 h)) H = 0, by (49), we obtain 



A Jcj) 2 h 2 dn + J g ij d p ,(4>h)d pJ (<j)h)dfi = \ h 2 g ij d p i 4> d pi cf> dfj, - \ h 2 <t>T4>dpL. (50) 
Setting / = he~ E { 2 ~ v / 2 we get 

'/' m d pl (4>h) dv = I gVdp* W) 8 p3 (0/)e~ 1 dp dx 



+ <j) 2 h 2 g l3 d p *Ed p] Edii 

+ \ jg lJ d P - (^f^dpjEQ-'dp dx 

and integrating by parts in the last term we get 
i jg i ^d p i{cj> 2 f)d pj EQ- 1 dpdx = ~ J^ 2 h 2 d p ,(g^d pl E)dfi>-oj J^h 2 d^ 

The identity (50) leads therefore to the inequality 

(A-w) Jtfh 2 d^i+J J <j?h 2 g ij d p >Ed p] Edix< j h 2 g 1 ^ 'fi^fl^d/i- / ' h 2 cj)T(j)d^. 

Let k = (ki,k 2 ) £ N 2 and 4> = 4> k {x,p) = 4>{x/ki)ip(p/k 2 ), where V € C£°, < 
ip < 1, tp = 1 on B(0, 1/2) and supp ip C B(0, 1). We obtain, denoting by C any 
positive constant, 

\h 2 d M + - j 4>lh 2 g l] d pl Ed p] Edn <^ f h 2 sup \g ij \x M< k 2 dn 

+ \{<f> k hV p E -V x (j> k ,h) H \ 
+ \{(j)khV x V- V P <f>k,h) H \. (51) 



(X — u) / ^ 
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Using Young's inequality, we can estimate the last two terms of (51) as 
| (<f> k h \7 p E • V^ fe , h) | < £ (-^- y^/i 2 |V^| 2 d/i + ej Jh 2 dp 
| {4> k hV x V ■ V p fc , h) | < f-L / ^d/x + e 2 / /i 2 d^ , 



for all ei,e 2 > 0, where C(^i) = su P|a;|</ i ; 1 {|Va : V'|}. Taking ei = C/{6k\) in the first 
line, e 2 = C£(A;i)/(4fc2) m the second line and using (47), we get 

(A-w-1) Ulh 2 dtx<^ fh 2 sup\g^\ XM<k2 d^ + C[^ + ^-) fh 2 dp. 

We see that ft = 0, taking first the limit fc 2 — )• oo and then fci — > oo. This concludes 
the proof of the theorem. □ 

B Finite propagation speed of relativistic kinetic 
equations 

This appendix is devoted to prove a general result that can be used to establish the 
finite propagation speed property for all relevant relativistic kinetic equations. It is 
obtained by adapting the proof of a celebrated uniqueness theorem for non-linear 
wave equations due to Fritz John [26], see also [34]. 

Lemma B.l. Let p,j G C 1 verify 

d t p + V ■ j = 0, t > 0, x e R d , (52) 

and \j\ < p. If p(0,x) = 0, for \x — xq\ < to, then p{t,x) = 0, for (t,x) <G A^q^q), 
where 

A(to, aso) = {(t, x) e [0, t Q ] x R d : \x - x Q \ < t - t}. 
Proof. Consider the function 

$(s,z) =t - [(t -s) 2 + to 2 (2t oS -s 2 )\x-x \ 2 ] 1/2 . 

Note that 

$(0,30 = 0, hm $(s,x) = t - \x-x \, and $|,._. oNto = 0. (53) 

Moreover, denoting R s (t , xo) = {(t, x) : t < $(s, x), |cc — a;o| < to}, we have 

A(to,-^o) = l>a<s<t Rs(t ,x ). 

Next we define 

p n (s,x) = p($(s,x),x), j n (s,x) = j($(s, x), x). 
Since p, j satisfy (52), then p n , j n verify 

d s p n = -V • j n 9 s $ + d s j n • Vx<t>- 
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Therefore, using (53), 

/ p n (s, x)dx = / d T p n (r, x) dr da- 

J\x-xo\<to J\x — x \<t JO 



I f (-V-j n a T $ + c» T jn-V K 0)drda; 

J \x— Xq I <£[} J 

- / V ■ j n <S>(s,x)dx 

J \x-x \<t 



whence 



\x— a 

in ■ W x (/)(s,x)dx, 

X — Xq I <t 



/ {p n -j n -\7 x (j))dx = Q. 

J \x — Xn\<tn 



' |a;-a;o|<to 

Moreover it is easy to verify that |Va;0(s,a;)| < 8(s) < 1. for all < s < to, thus, 
since in addition \j\ < p, we get 



Pn da; = =>■ p = on A(to, a;o)- 

X — Xq I <t 



□ 



The preceding lemma can be applied to any relativistic kinetic equation which 
is compatible with the continuity equation (52). Precisely, to any kinetic equation 
of the form 

d t f+p-V x f = Q[f], 
where Q is a (possibly non-linear) operator such that 

Q[f}(t,x,p)dp = 0. 

i 

The previous identity implies that 

/dp, j = fpdp 



satisfy the continuity equation (52). Moreover, since \p\ < 1, then |j| < p and 
Lemma B.l applies. 
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